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Abstract 



We present a general formalism for computing Cosmic Background Radi- 
ation (CBR) and density fluctuations in open models with stiff sources. We 
decompose both the metric fluctuations and the fluctuations in the stress- 
energy tensor into scalar, vector and tensor modes. We find analytic Green's 
functions for the linearized Einstein equations in the presence of stiff sources 
and use this formalism to estimate the amplitude and harmonic spectrum 
of microwave background fluctuations produced by topological defects in an 
open universe. Unlike inflationary models that predict a flat universe and a 
spectrum of CBR fluctuations that is enhanced at large angular scales, defect 
models predict that CBR fluctuations are suppressed on angular scales larger 
than that subtended by the curvature scale. In an = 0.2 — 0.4 universe, 
these models, when normalized to the amplitude of CBR fluctuations observed 
by COBE, require a moderate bias factor, 2 — 3, to be compatible with the 
observed fluctuations in galaxy counts. In these models, accurate predictions 
can be made which are testable through CBR experiments in the near future. 
A CBR measurement of $7 would then be possible, up to the limit imposed 
by cosmic variance. We discuss some of the philosophical implications of an 
open model and propose a solution to the flatness problem. 



I. INTRODUCTION 



In recent years, most theoretical work in cosmology has assumed that the universe is flat 
and matter-dominated. This assumption is nearly inevitable in inflationary scenarios that 
predict that Q should be close to unity. However, defect models make no such predictions 
for the density of the universe. 

There is a host of astronomical evidence that suggests that the universe may be open. 
White et al. |^ argues that X-ray observations of clusters imply that at least 20% of their 
mass is baryonic. When combined with standard hot big-bang estimates of the baryon 
density [^, Qbh"^ ^ 0.015, this implies that the total density in non-relativistic matter is 
much less than unity. Observations of galaxy random velocities [H flnd that as — 300 km/s, 
a factor of three below the predictions of COBE normalized flat scale-invariant cosmologies. 
The matter density inferred from comparisons between the galaxy correlation function in 
redshift and real space is also much smaller than unity and is compatible with Q ~ 0.3 [Q. 

Over the past decade, observational cosmologists have devoted much effort to measuring 
the spectrum of density fluctuations and have found signiflcant evidence for more large- 
scale structure than predicted in flat universe models. As the predicted spectrum of density 
fluctuations is peaked on the physical scale corresponding to the horizon size at matter- 
radiation equality, ?,2/VLh? Mpc, models with f2 = 1 predict less large scale power than 
low Vt models. For example, the standard CDM model underpredicts the ratio of galaxy 
fluctuations observed on the 30/h Mpc scale in the 1.2 Jy survey to galaxy fluctuations 
inferred on the 1/h Mpc scale in the same survey by nearly a factor of 5! On the other hand, 
CDM models with VLh ~ 0.2 are remarkably successful at fltting observations of large-scale 
structure f^. 

In a flat universe, topological defect models also fail to produce the observed large- 
scale structure. Albrecht and Stebbins found that the predicted spectrum of density 
fluctuations in a cosmic string model is peaked on the string coherence scale, which is even 
smaller than the horizon size. Because of this lack of large-scale power, they concluded that 
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a string-seeded CDM model was not compatible with the observed large-scale structure and 
failed more dramatically than the inflationary scenarios. In Pen, Spergel and Turok [|^, 
hereafter PST, we explored global defects such as global strings, monopoles, texture and 
non-topological textures. While monopoles, texture and non-topological textures all predict 
a power spectrum of density fluctuations with more large-scale power than the cosmic string 
model, these models still fail to produce the observed level of galaxy fluctuations by nearly 
an order of magnitude on the 20/h Mpc scale. This failure is due to the defect coherence 
scale at equality being too small. As this scale is proportional to flh"^, this problem will be 
alleviated in a low Q universe. 

In the excitement following the announcement p,^] that COBE has detected thermal 
fluctuations in the cosmic microwave background (CBR), many cosmologists declared that 
the detection was evidence for a flat universe and the inflationary scenario. This, however. 



need not be the only interpretation of the COBE results. The two-year COBE results [21 



are 



also consistent with open universe models with either adiabatic density fluctuations |10,11 



or equation-of-state density fluctuations [in 



Inflationary scenarios predict either a scale-invariant spectrum of CBR fluctuations or 
a spectrum of fluctuations that has more power on large angular scales []T^. In models 



with signiflcant gravity wave contribution [|I^ and in power law inflationary models [|18 



the low multipole fluctuations are enhanced relative to fluctuations on small angular scales. 
Inflationary models with cosmological constant [|l^] also predict enhanced contributions on 



large angular scales. On the other hand, analysis of COBE's two year data set p0|j21 
suggest that the low multipoles are not enhanced but suppressed relative to the fluctuations 
on smaller angular scales. While there is still significant statistical uncertainty in the two- 
year COBE data, the Tenerife CBR observations [|I^ and the FIRS results [0 also hint that 
the slope of CBR fluctuations may be steeper than predicted in any inflationary models. 



This steep CBR fluctuations spectrum is one of the predictions of open universe models []T^ 
If the suppression of the low multipoles is still seen in the four-year data, then the COBE 
data may turn out to be incompatible with most inflationary scenarios. 



If the universe is open, then topological defect models are particularly attractive. The 
basic concept of inflation lies in solving the monopole and horizon problems through an 
extended de-Sitter phase. As other authors [|T5| have pointed out, the curvature of space- 
like section in a de-Sitter and also an empty universe depends on the choice of coordinate 
systems, and could be either fiat or hyperbolic. The prime reason for considering a flat 
universe a generic prediction of inflation is the presence of the flatness problem. In this 
paper we will propose a solution along the lines of Linde, by using the weak anthropic 
principle as a selection effect. It hinges on the huge photon to baryon ratio (approximately 
10^'^), and gives a natural scale to the problem. The weakest point of inflation has been the 
fine tuning of parameters required to generate the observed potential fluctuations. So even 
within the frame work of inflation, topological defects may still be desirable as they explain 
this fine tuning more naturally. 

But quite independently of the existence of inflation, if we are left to find another route 
for explaining the large-scale homogeneity of the universe, it would be most natural to 
assume that the universe started with smooth initial conditions and that causal physics 
generated density fluctuations. Spergel |T^ presents analytical arguments that suggested 



that defect models are likely to be compatible with observations of large-scale structure in 
an open universe. In this paper, our numerical calculations support these arguments and 
suggest that these models merit careful consideration. 

This article builds upon the analytical and numerical techniques presented in PST. In 
PST, we attempted to quantitatively compare the predictions to topological defect models 
in a fiat universe to the COBE observations and observations of large-scale structure. This 
paper extends our approach to an open universe. Our GBR calculations make no assumptions 
about the nature of the dark matter. However, when we compare the COBE normalized 
theories to the observed large-scale structure, we assume that the universe is dominated by 
cold dark matter. 

In this paper, we estimate the amplitude of GBR fluctuations generated by various 
topological defect scenarios. Because of ease of computation, we focus on the texture models. 



While we have deferred the challenge of evolving a string network in a hyperbolic universe, 
we believe that qualitatively the results of this paper can be extrapolated to other defect 
models. In section [I^, we extend the analytical formalism that we developed in PST to open 
universe models. In appendix |B|, we also extend the formalism to vacuum-dominated models. 
In appendix |C| we extend to formalism to closed models. In section [m| , we describe our 



numerical algorithms for computing defect evolution. In section we present our numerical 
results and emphasize the characteristic CBR signature of defects in open universe models 
and then we discuss the predictions of the COBE^normalized open universe models for large- 
scale structure. In section ^ we discuss philosophical motivations for open universes. In 



section VI, we sum up. 



II. CBR FLUCTUATIONS IN AN OPEN UNIVERSE 

In this section, we modify the formalism developed in PST so that we can calculate the 
amplitude of CBR fluctuations produced in an open universe. 

The assumptions entering in the calculations are the weak field limit for gravity, which 
allows us to treat gravitational perturbations in linearized form. The quantum field is 
evolved as a classical field, as it is a boson field with large occupation number. An additional 
assumption is the stiffness of the source term, that gravity does not affect the evolution of 
the defects. This is certainly justified for global defects, which dissipate their energy into 
Goldstone modes. In the case of gauged cosmic strings one needs to account for their 
energy dissipation in gravity waves. The defect field is governed by the nonlinear a model, 
which is a highly nonlinear evolution equation. While certain scaling laws can be computed 



in a fiat space time |T6[, this is not possible in the transition regime between matter and 
curvature domination. The field correlation length can not be accurately specified, and most 
analytic approaches are no longer applicable. Thus one needs to simulate the field evolution 
numerically, and measure these quantities. 

The gravity is, however, still linear, and we will show below how to solve that problem in 
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an open universe once we are given the source terms. They are still straightforward integrals 
over the energy momentum tensor with certain Green's functions. 

In PST, we decomposed variations in the source stress energy tensor and fluctuations in 
the metric into scalar, vector and tensor components. Variations in the trace of the spatial 
stress energy tensor generates growing density modes that can form galaxies. The traceless 
scalar source term, the vector source term and the tensor source terms do not generate 
growing density modes, however, they do source decaying metric fluctuations that produce 
CBR fluctuations. We have reexamined our flat universe simulations and found that the 
CBR fluctuations generated by the scalar modes are the dominant source of fluctuations: 
the scalar growing mode term alone accounts for 70% of the CBR fluctuations. 

The decomposition into scalar, vector and tensor modes is a non-local calculation that 
is numerically challenging in an open universe. In order to evaluate the viability of defect 
models in an open universe, we will focus on only the contributions of the growing scalar 
mode in this calculation. In an open universe, we expect that the vector and tensor modes 
are even less important than in a flat universe as they are suppressed relative to the scalar 
modes by powers of (v/c), where v is the defect velocity. In an open universe, the rapid 
expansion of the universe slows the defect velocities. Because of our ignoring the anisotropic 
stress, vector and tensor modes, the amplitude of CBR fluctuations calculated in this paper 
should be multiphed by a factor between 1 — 1.2. 

In PST, we showed that the variations in the trace of the spatial stress, ©, sources 
variations in the scalar piece of the metric [PST 46] : 



h- + 2-h- = -STfCe = S 



(1) 



a 



These metric fluctuations contribute to the Sachs- Wolfe integral: 




(2) 



where 




h 



(3) 
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[PST 50] and 



j + 2^-J = -\h- 
a 6 



(4) 



[PST 49]. Here, a is expansion factor and dot denotes derivative with respect to conformal 
time. In the notation of PST, J = k'^J'. 

After a httle algebra, equation and can be combined to yield an equation for 

the scalar stress contribution to CBR fluctuations: 



K{x,ri) = / driS{x,ri) 
Jo 



a{riY d{ri) ,^.2r-/ ~\ 
ayriY a[ri) 
( a{ri) dijj) \ aijiY 



_ 3 

a{vi) (i{riy I a{r)y 



(5) 



where 



V drj 



a{r]y 



(6) 



and 



H{r],r]) = / a{r]) dr] / 



n dr( 



(7) 



aivCY 

Note that E can also be used to relate the vector stress energy source term to the vector 
metric fluctuations term [PST 47] and to evolve the decaying scalar metric fluctuations [PST 
49]. 

In a flat matter-dominated universe, aij]) = rf, thus, E{r],fj) = l/3r]^ — l/3fj^ and 
if(?7,f/) = 7]^ /Ibff — rf /Q + rf /IQ. Combining these results with equation yields: 



K{x,ri) = dfiS{x,r]) \^ 



(8) 



This simple result is due to the simple form of a{7]). 

In an open matter-dominated universe, the expansion factor has a more complicated 
form: 



a{ri) 



-/C 



cosh(V— /C??) — 1 



(9) 



where flo is the density in matter today and JC is the curvature scale. For the rest of the 
paper, we set /C = —1 and use it as the physical length scale in the calculation. Note that for 
small 1], equation approaches the fiat space form. Combining equation (|^) with equations 
(I) and (0) yields: 



and 



TT/ _s 1 5 

a b b 6 



a d{ri) [1 — d] 



(10) 



bd 



n2 



a{r]){a - 3) + 8(77 



(11) 



where d = a{fj) and a = a{r]). Equations @, (|T0|) and ([TTD can now be combined to yield 



rrj _ i (l^ 

K{x,ri) = / dfjS{x,ri)< — 
Jo a'^ 



^ 3 ^U^ /\¥~3~6, 



+a(l — a) 



+ f ^ + 2^ '^(^ - 3) + 3r7 
3a \a J 6a^ 



+ 



a — a(l — a)ri 



5 + 2a 



(12) 



In the limit of small rj, equation ([T^) reduces to equation (^. 

The vector modes still allow a simple integral, for which we have 



hY = ^ I aWf^fdv'. 



(13) 



Tensor modes, on the other hand, propagate and require two Green's functions to express. 
As in [PST 54], we write 



hJi = IbnG 

''J 



where now the Green's functions are 



G2 



W{ri') 



cos{kr]) sm(kr])d 

a 2a^ 
cos{kr])d ^ sm{kr]) 



df)' 



2d^ a 
W = G1G2 - G2G1 



(14) 



(15) 
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We have assumed a decomposition in terms of eigenf unctions of the Laplace-Beltrami op- 
erator V^, which in our framework will be the sum of two sine frequencies and one Bessel 
function index. Note that the tensor mode Green's function is equation (^) are valid for 
open, flat and closed models. 

In the next section, we present our numerical algorithms for evolving the defect field to 
compute S{x,'fj), of integrating equation (|12D to compute the metric fluctuations and for 
following photon trajectories to integrate equation (||). 

III. NUMERICAL IMPLEMENTATION 

In this section we describe the numerical implementation issues. The orig- 
inal program is freely available by anonymous ftp from astro.princeton.edu in 
/upen/StiffSources/openuniverse. It is written in standard C++ and C, and should compile 
and execute on any machine with these compilers. It is optimized to execute very efficiently 
on the convex vector architecture. In fact, simulations are always limited by memory, be- 
cause the volume of a hyperbolic universe is exponentially large, so the calculation scales as 
0(A^log(A^)) where is the memory requirement. 

The basic strategy will be to apply the mode decomposition from PST. In order to work 
in an open universe, many changes need to be applied which are described below. 

A. Grid 

The very first obstacle is the formulation of a regular lattice to discretize a hyperbolic 
manifold. The requirements are: 1. it must have constant volume per lattice element, 2. 
appear locally Euclidean, 3. be easily mapped onto the serial storage of a computer, 4. allow 
the Laplacian to be easily invertible. 

For this purpose the Poincare metric provides a very nice tiling, which retains many 
of the regularities of a fiat space Cartesian lattice. The spatial metric is given by the line 
element 
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2 dx^ + dy^ + dw^ 
ds^ = . (16) 

With a change of variables z = ln{w), we satisfy all the requirements stated above. On 
small scales, it is explicitly Euclidean, which simplifies the implementation of the differential 
operators. This metric maps onto the more familiar Friedman coordinates through the 
change of variables 

cos(0) sinh(x) 

cosh(r7 — x) 
sin(0) sinh(x) 
^ cosh(?7 — x) 
cosh(r7) 
cosh(r7 — x) 

cos{9) = tanh(r7) (17) 

The mapping onto the Friedman coordinates is shown in figure ^ The inverse mapping 
is displayed in figure The salient features of this metric are the explicit translational 
symmetry along all three dimensions, and the rotational symmetry about the z axis. The 
only explicit numeric anisotropy occurs for rotations in the z — x and z ~ y planes. This 
effect is easily tested for in the simulations by checking the alignment of the quadrupole 
with the coordinate grid. 



B. Technical Issues 

For simplicity, we will work in units where At = Az = 1. The horizontal discretizations 
exp{—z)Ax, exp{—z)Ay are adjusted to be as close to unity as possible, while still satisfying 
the periodicity constraints. In these units our free discretization parameters are the curvature 
radius R, the mesh height H and the periodicity length L at z = 0. We choose —H/2 < 
z < H/2. We divide the computational grids into tiers Tj at constant Zi, each of which 
is a square matrix. We need to subdivide into an integral number of lattice points, for 
which we calculate the integer such that the area bounded by each lattice point most closely 
approximates dz"^. Since the Tj are represented by regular matrices, all parallelizations and 
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vectorizations are performed at this level. The tier concept is then implemented as a C++ 
class. The key operation is the projection (interpolation) of one tier into the geometry of 
its neighboring tier above or below, 

P+ : T, T,+i, P_ : T, ^ T,_i. (18) 

In order for the discretization errors to be small, one needs to have several grid cells per 
curvature length. Another limiting constraint is the scaling behavior. In flat space, many 
defects achieve a scaling solution where the energy density is proportional to 1/a^. We 
require the numerical solution to achieve such a scaling law before the horizon size grows to 
the curvature scale. 



C. Mode decomposition 



Our basic tool is the Fast Fourier Transform (FFT). Since the grid is periodic in x and 
y, we can write any function as a sum of Fourier components, 

2i7[{nx + my) 



^exp 



L 



(n, m, z). 



(19) 



In order to retain discrete orthogonality, the numerical grid points must be all aligned at 
the phase origin x = y = 0. Then we simply keep a different number of Fourier modes at 
each z. The Laplacian then becomes a second order ordinary differential equation in z, 



An'^in^ + IV?) 



L2 

which we integrate to second order accuracy. 
In terms of the discretized variables, 

cos{7Tn/nk) + cos(7rm/mfc) — 2 



dz dz' 



(20) 



+ 



e-x.^{-Zk) 



exp(2;fc+i/2)(0fc+l - 0A:) - exp(^fc+i/2)((/)fc - 



(21) 



Note that the maximal mode n, m depends on the level of the tier. In the case that 
the tier above does not contain a corresponding mode, we use a zero boundary condition. 
Equation (^Tj) is a tridiagonal system which is solved in linear time. 
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Some care needs to be taken with the boundary conditions. One can easily violate 
causality from the non-local inversion of the Poisson operator. We thus need a boundary 
condition consistent with causality. A simple approach would be to set the boundary at the 
edges of the computational domain, Zu, zi to zero. 

In analogy with electro-magnetism, the zero boundary conditions can be physically in- 
terpreted as a distribution of surface charges which cancel the desired fields. They cause 
waves to be reflected at the boundaries, so that the gravitational field boundary conditions 
become consistent with the field evolution. 

The moving boundary condition can source scalar field and gravity waves, but since these 
only travel at the speed of light, a buffer zone will prevent them from affecting the photon 
cone in the calculations. 

D. Field Evolution 

We implement a nonhnear a model following the same approach as used in PST. The 
equation of motion for the continuum field reads: 

\dyd^(l) = + A0 (22) 

where A is a Lagrange multiplier which must be chosen to satisfy the constraint (p^ — 1. 

The time discretization has two degrees of freedom, corresponding to the initial values of 

and 0, which we represent through the field configuration at two consecutive time steps. 

We proceed in two steps. First we calculate the Laplacian. Then we advance the field in 

the direction of the Laplacian subject to the constraints of the nonhnear a model. We treat 

the two issues in turn. 

In the Poincare metric, the Laplacian is expressed as 

2^' z z 

= exp(-— )(0,^^ + (p^yy) + exp(-— )a^exp(— )0,^. (23) 

The first two terms are trivial to calculate with the standard central difference formula. 
With help of the projection operators P+- we can easily evaluate the vertical derivative 
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exp(--^) 



exp 



(24) 



Along the x, y axes we can simply use periodic boundary conditions. At the top of the grid 
z — H/2 we simply extend our grid upward, which costs very little in computational effort 
or memory because very little volume is enclosed in that region. The bottom boundary 
needs to be treated more carefully. We choose a safety buffer zone of a few grid cells below 
the last grid point that is traversed by photons. 



E. Implementation 

For texture models, there is reason to believe that the primary contribution comes from 
the spatial trace of the energy momentum tensor, 

e^jZTl (25) 
1=1 

In flat and empty space, the exact texture solution generates an energy momentum tensor 
which is such a pure trace. As we describe below, the problem is isomorphic to such a flat 
problem in both the early and late time limits. The main contribution thus arises from the 
trace. In our current implementation, we have chosen the approximation to only retain the 
trace part. For other defects, especially the cosmic strings, other components of the energy 
momentum tensor are expected to play a dominant role and one needs to implement the full 
mode decomposition described above. 

The simulations are constrained on several ends. We need a fair number of grid cells 
per curvature radius in order to achieve an accurate flat scaling density for the fleld before 
the numerical horizon size becomes comparable with the curvature scale. A violation of this 
constraint would cause the simulation to enter the curvature transition with an incorrect 
energy density, which would appear as a systematic subsequent error. In practice this 
corresponds to about 8 grid cells per curvature radius. 

With the 2 Gigabytes of memory on our convex C3440 we can run simulations down to 
O = 0.2, which take about 4 hours of (wall clock) execution time. 
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F. Tests of the Code 



The two extreme limits of the parameter space have exact solutions. When 17 = 1, we 
have an expanding fiat space and we can test the exact scaling solution for a single unwinding 
texture as we did in the fiat space calculations [0. The other limit f2 = is an empty 
universe, which is nothing more then Minkowski space. Using the Milne transformation, we 
map the exact scaling solution for a texture as initial condition on our grid, and test for the 
subsequent evolution. In an empty universe using conformal coordinates, the cosmic scale 
factor a = 1], so the Milne solution to the Einstein equation in an empty universe becomes 

rfs2 ^ _^^2 _^ ^2^J^ ^ ^2^^2) (26) 

1 + q'' 

With a change of coordinates r = qr, t = ry/1 + t^^q recover the Minkowski metric 
ds"^ = —dt^ + c/r^ + r'^dVl^. Our numerical grid performs very well on this test. Since the 
same code has very small errors on both limits, we can claim some confidence that it should 
perform well in between. 

IV. RESULTS 

A. CBR fluctuations 

Using the algorithm outlined in the previous section, we have computed the CBR fluc- 
tuations produced by scalar potential fluctuations in an open universe. The results of these 
calculations are shown in figures |^, ^ and |^ for four different observers in 17 = 0.21,0.4, 
and f2 = 1 universes. The anisotropy of the grid is almost visible as an enhancement of the 
quadrupole in figure ^for Vl = 1. In this case the light rays always move at a constant angle 
to the grid, so we conclude that even in the worst case, the grid anisotropy only has a minor 
effect. 

In an open universe, the defects dynamics slows down due to the rapid expansion of the 
universe, which exhibits itself as a loss of power on large angular scales, in particular the 
quadrupole terms. 
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The different lines in the figure denote the results from different realizations and the 
spread in values is a measure of the variance in c;. The abscissa in the plot is the amplitude 
of the multipole moments, q = Z^m^^Li/l^^ + 1), weighted by l{l + 1). In a fiat inflationary 
model with scale-invariant spectrum of fluctuations, q/(Z + 1) is a constant for / << 200. 
Our results imply that for defect models in a low Q universe, the shape of the multipole 
spectrum is qualitatively different from other models. The low order multipoles are strongly 
suppressed, while the multipoles on scales significantly below the angular scale subtended 
by the curvature scale today still have a scale invariant form. 

For purposes of comparison with the DMR measurements of the temperature fluctuations, 
we have fit the results of our numerical simulations with a fitting form: 



Cll{l + l) = Co/[l + {lmax/l) 

The simulations are fit by Imax — Sfl^'^. Even if the universe is flat, the quadrupole is 
somewhat suppressed in any model with topological defects as defects on scales comparable 
to and greater than the horizon size have not yet had time to collapse. 

The qualitative features of topological defects in an open universe can be approximated 
using any flat universe calculation. To first approximation, one can consider the photon 
sphere we observe today projected back at a redshift of 1 + z ~ 1/Q. This will also produce 
a flat spectrum for large I, with a white noise cutoff for small I. We used our flat space 
stopped at z = 1 and z = 1.5, and the multipole spectrum is depicted in figures ^ and |^. 

The COBE two year observations have only been analyzed under the assumption of a 
power law spectrum of CBR fluctuations [^,^,^. These analyses conclude that COBE 
has measured a quadrupole of only 6 ±3/iK while their fit to larger multipoles imply Qrms = 
20/iK. In an inflationary model, such a small value for the quadrupole should be observed 
in less than 5% of the universe. In a defect model in an open universe, the low value for the 
quadrupole is predicted. On the other hand, COBE two-year measurements do not find a 
suppression of the / = 3 and / = 4 modes. 

As discussed in section 3, our calculations only include the contributions from the scalar 
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growing mode to the CBR fluctuations. The calculations do not include contributions from 
decaying modes, vector fluctuations and gravity wave fluctuations. These contributions are 
sub dominant in a flat universe and should be even smaller in an open universe. If these 
terms were included, then the amplitude of the CBR fluctuations would be increased by a 
factor / ~ 1.0 — 1.4. The upper bound is based on our flat space calculations [0. Thus, the 
amplitude of CBR fluctuations needs to be multiplied by this factor. 

Topological defect theories have one free parameter: the scale of symmetry breaking, 
00- Note that the abscissas in flgures H - H need to be multiplied by Stt^G^q to convert 
the results of the calculations into temperature fluctuations. We will flx this parameter by 
normalizing our results to the COBE two-year observations by convolving our results with a 
Gaussian beam with full width half maximum of 10° and flxing {ST/T)f.^^ to the value of 40 



/iK suggested by harmonic analysis of the two year data . The second column in table 1 
shows the normalized value of Stt^G^q for different values of Q. Note that these values are 
slightly higher than our earlier calculation based on the DMR one-year data. 

While our calculation did not compute the amplitude of CBR fluctuations on small 
angular size, we can extrapolate calculation on defects in a flat universe to the open case. 
Coulson et al found that in a re-ionized universe, the multipole spectrum was flat from 



large angular scales to / ~ 60. This multipole moment corresponds to the angular size 
subtended by a texture collapsing near the surface of last scatter in a re-ionized universe. In 



an open universe, the relationship between horizon size and angle is altered: 6 ~ 



5^-0.5 



for z » Q^^. Thus, we expect that in an open re-ionized universe, the multipole spectrum 
would be flat from / ~ Imax to / ~ 60Q~^'^. 

If the early universe was not re-ionized by a generation of star formation before z ~ 50, 
then CBR observations on small angular sizes are probing the universe at z ~ 1300, the 
epoch of recombination. In both defect models and scale-invariant curvature models, there 
should be a "Doppler peak" at / ~ 200^2^*^'^. In an inflationary model with curvature 
fluctuations, this peak is produced by the sum of velocity perturbations and the product 
of potential fluctuations with entropy fluctuations. In these models, the fluctuations are all 
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produced by the growing modes. In a defect model, the "Doppler peak" is due primarily 
to the entropy and potential fluctuations produced by the decaying modes excited by the 
collapse of defects [^. As in curvature models, this peak should occur at / ~ 200fi~°'^, 
thus measurements of its location may provide a determination of Q. 

Just as in a fiat universe, one of the distinctive predictions of a defect model are the 
non-Gaussian character of the temperature fluctuations. This non-Gaussianeity should be 
most apparent not in the distribution of temperature fluctuations, but in the distributions of 
temperature gradients [p4| , p5| . In a flat universe, Coulson et al. find that this non-Gaussian 
character is most apparent on angular size of ~ 3°, the angular size subtended by the surface 
of last scatter in a re-ionized universe. In an open universe, this angular size is shifted to 
~ 3i7°'^ degrees. 

While simulations were done only for texture models, the results can easily be extrapo- 
lated to other defect models in an open universe. In all of these models, we expect a similar 
suppression of the low multipoles as the rapid expansion of the universe slows the evolution 
of the defects responsible for generating fluctuations on large angular scales. 



B. Large Scale Structure 

Having normalized the defect model to fluctuations in the GBR, we now turn to predic- 
tions for mass fluctuations on the scale probed by galaxy surveys. Fluctuations on scales 
smaller than ~ 1000 Mpc were generated when Q was close to unity, thus, the results of our 
earlier work on density fluctuations in a flat universe can be directly extrapolated to open 
models. 

In Pen et al. [0, density fluctuation computed by numerical simulations flt by a function 
[PST 32, 33]: 

P(k) = ^(^)' 

Qh^ [1 + {/3k) + i-fky-^]^ 

where D{Q) ~ is ratio of the linear growth to that in anQ = 1 universe, a = 225(e/3.7x 
lO-'^Y/{nh^), (3 = 3.5/{nh'^) Mpc, and 7 = 2.75/{nh'^) Mpc. Using the normalization in 
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Table 1, this directly yields the amplitude of mass fluctuations predicted by linear theory. 

On the scale of tens of kiloparsecs, light clearly does not trace mass. This is the source of 
the missing mass problem. It is less certain whether light traces mass on the scale of several 
megaparsecs. Cosmologists parameterize this uncertainty by a bias parameter, b, the ratio of 
the variance in the fluctuation in the galaxy counts to the fluctuations in the cosmic density 
fleld. Here, we determine the value of b needed to flt astronomical observations. 

The QDOT survey measured the fluctuations in galaxy counts by obtaining redshifts to 
a large infrared selected galaxy sample. Saunders et al. smoothed their galaxy density 
fleld with a Gaussian smoothing window with fllter length of 5h^^, lOh^^, and 20h^^ Mpc 
and found a variance of 0.436 ± 0.091, 0.184 ± 0.05 and 0.0669 ± 0.019 in the density fleld. 
The fourth, flfth and sixth column in Table 1 lists the required bias factor needed to flt the 
central value in the QDOT survey. The statistical uncertainties in the QDOT survey and 
the COBE measurements lead to a ~ 25% la uncertainty in b. This is in addition to the 
uncertainty due to the limitations of our numerical calculations. 

In flgure ^ we compare the predicted power spectrum of density fluctuations to the power 
spectrum of galaxy fluctuations inferred from APM survey. In this flgure, we assumed a bias 
factor of 1.5. 



V. SPECULATIONS 

At various points in history, different choices for the curvature of the universe have 
been considered most natural. Einstein initially considered an eternal and flat universe 
with a cosmological constant most aesthetically pleasing. But with Bubble's discovery of 
the universal expansion, the common belief was that the universe should be a closed three 
sphere, which is bounded in both space and time. This universe would end in a few Hubble 
times, and again we would live in a very ordinary epoch. We will use the same Copernican 
principle to argue that an open universe is almost as well suited. 

In the last decade it has become fashionable to return to consider spatially flat universes 
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at most appealing, since there would not be any curvature scale which needs to be explained. 
This is primarily due to Dicke's anthropic argument, who used the Copernican principle to 
argue that we should not be living just at he end of the flat epoch. Not to forget that the 
inflationary paradigm, which appeals to the de-Sitter model to solve the so-called horizon 
problem, is simplest to explain in a scale-free scenario of a flat universe. But the very 
problem that they are invented to solve, the absence of a preferred curvature scale, leads 
to a clash with the Copernican principle as we live at a very special time, just near the 
beginning of a matter-dominated universe which would now last for a truly lengthy period 
of proper time. As we have seen, in the absence of a perfect fluid to label a preferred 
coordinate frame, the curvature of a universe can be transformed from flat to hyperbolic by 
a gauge choice. The same holds true for a de- Sitter space. 

An open universe model does have aesthetic advantages that have been at times over- 
looked. In an open universe, we are most likely to live in the brief period of time between 
radiation and curvature domination. Density fluctuations do not collapse during the radia- 
tion dominated epoch and are growing logarithmically slowly during curvature domination. 
If a universe went directly from a radiation dominated phase into curvature domination, no 
structure or life would ever be conceivable. It is only due to a lucky coincidence that we 
have a slight baryon asymmetry of 77 = 10~^^, possibly due to baryogenesis in the electro- 
weak phase transition. This allowed structures to form through gravitational instability in 
the short interval between matter-radiation equality and curvature domination. Since this 
interval lasts only for a finite time, and we live approximately in the center of this period, 
there is no violation of the Copernican principle. We thus appeal to the observed smallness 
of electro-weak baryogenesis to set the scale for the hyperbolic curvature. Dirac's small num- 
ber is not really a single small number. The radius of curvature and the horizon size today 
are simply the product of the proton mass, the baryogenesis photon to baryon ratio, and 
the smallness of initial fluctuations, observed by COBE to be 10~^, which in the topological 
defect framework arises from the GUT scale. Such a model has no fine tuned parameters. 

From a geometric viewpoint, a closed universe is appealing due to its simplicity: a 
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three sphere is the unique universal covering of positively curved three-manifolds, and has 
a finite volume which would certainly be an attractive property for any designer or process 
which might have created the universe. But it is interesting to note that there are only a 
finite number of alternate global topologies which such a designer has to choose from. The 
projective 3-sphere is one such example. If one analyzes negatively curved spaces, one 
can of course consider the global covering H^, but there are many alternatives, including 
the periodic Poincare space which we utilized in this paper. The name "open" only applies 
to the local properties, and we can certainly have a spatially closed "open" universe. If one 
considers only hyperbolic spaces of finite volume, one finds an infinite number of possible 
topologies. At fixed curvature R^, the volumes of these topologies can have collection points 
on the real line, and one might expect our universe to be chosen from a topology near such 
a collection point. A number of authors have attempted to calculate transition probabilities 
between these configurations in 2+1 dimensions These studies suggest that topological 
change may well be possible. Unfortunately, as with most quantum gravity calculations, 
many infinite quantities arise in the process, making is difficult to uniquely predict the 
outcome of such an estimate. 

Whether the negative curvature results from quantum gravitational tunnelling, or an 
alternate exit from the de Sitter phase, or some other yet unknown means, the anthropic 
principle does set a minimum scale to the curvature radius. In order to bound it from above, 
one could argue that the intrinsic process forms hyperbolic space-times with small curvature, 
most of which are not observable. So we might live in the smallest allowed scenarios which 
allow nonlinear structures to form. 

However, as we lack a theory of quantum gravity that can predict whether a flat, open or 
closed universe is most likely, we believe that all of these models merit careful consideration. 
Ultimately, this question must be resolved observationally. We have shown that the philo- 
sophical arguments that are been invoked to argue for a flat universe are quite ambiguous, 
and could just as well be used to argue for an open scenario. 
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VI. CONCLUSIONS 



In this paper, we explored the cosmic microwave background signature of the formation 
of large scale structure by defects in an open universe. We have described a new efficient 
exact solution of the linearized Einstein equations in hyperbolic FRW spacetimes. The 
otherwise expensive mode decomposition can be implemented very efficiently thanks to 
the Fast Fourier Transform in the Poincare metric. This formulation was then applied to 
calculate predictions of texture models in an open universe. We then addressed the classical 
philosophical arguments including the flatness problem. We showed that the same anthropic 
and Copernican arguments that were used to argue for a fiat universe, are in fact better 
satisfied in an open model. The curvature scale is naturally explained as a product of three 
moderately big numbers, which are experimentally well established. 

For Q = 0.4 and h = 0.7, the power spectrum of density fluctuations in a COBE- 
normalized texture model has the correct spectral shape and is consistent with the observed 
level of galaxy fluctuations for 6 = 2 — 4. The uncertainty in normalization is a combi- 
nation of the numerical uncertainties in our calculations and the statistical uncertainties 
in the observations. A model with fl = 0.4 and 6 = 2 is consistent with various dynam- 
ical measurements of Q on the scale of clusters and superclusters. The observations that 
are most problematic for this model are the large-scale streaming velocities inferred from 
various proper distance surveys. 

While our work focused on textures in an open universe, we expect qualitatively similar 
results for other defect models. The basic results appear to be governed by geometry and the 
changed relationship between angle and physical scale. This is apparent when we compare 
a flat universe simulation stopped at 1 + z — and rescaled by a factor of Q in angle with 
an open universe simulation. 

Defects in an open universe make a distinctive prediction for the CBR spectrum. In 
these models, very few fluctuations are generated at late times and at large angular scales. 
Thus, the models predict a suppression of the quadrupole and other low multipole moments. 
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The low value of the quadrupole detected by COBE is consistent with this set of theories. 
A detailed analysis of the COBE DMR results is needed to determine which range of val- 
ues of fl are compatible with the observed universe. A qualitative understanding is quite 
straightforward and is obtained by rescaling the flat space spectra and introducing a break 
near the curvature scale. 

Defects in an open universe are a viable alternative to popular scenarios for structure 
formation and merit closer study. 
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Geometric Formulae in Hyperbolic Spaces 

APPENDIX A: NOTATION 



Rc will denote the curvature radius of the universe, which is given by Rc = c/ {Ho^/l — Q) . 
A subscript of denotes a parameter's present value. We use t to denote the proper time, 
7] to measure conformal time. They are related through 

t = — ^ — {smhrj — T]) (Al) 
sinh rjo 

In conformally hyperbolic coordinates, the expansion factor is a = coshr] — 1 = (2/^2) — 2. 
The conformal time measures the number of comoving curvature radii traversed by a photon. 
It is given by rj = cosh"^(2/r2 — 1). We write the FRW metric as 

ds^ = -dt^ + aitfi^^ + r^dn^). (A2) 

Proper distance between two time synchronized observers is As = ai?c sinh~^r. The coor- 
dinate r is in units of comoving curvature radii. 
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APPENDIX B: VACUUM DOMINATED UNIVERSES 



This appendix contains results for a vacuum-dominated flat universe containing dust 
and vacuum energy. In this universe, the expansion factor can be computed from the energy 
equation: 



'da" 
'dt 



a aQ 



where f^o, clq and Ho are the density in matter, the expansion factor and the Hubble constant 
today and t denotes physical time. This equation can be solved to yield: 

a{t) = ^Q'J\l-Qoy^'smh'/\Hot) 

The conformal time, rj in a. flat vacuum dominated universe can be computed with the 
aid of Gradshteyn and Rhyzik equation (3.166.22): 

f"* dt r da 



V 



dx 



a{t) Jo a{t){da/dt) Jo + x^) 

_^ f 1 + {1 - VS)a\ \/2 + V3 



3V4 



cos 



(Bl) 



,l + (l + x/3)ay ' 2 

Here, is the elliptic integral of the flrst kind. 

It would appear that this complicated relationship between a and t] would make it im- 
possible to evaluate equations (||) and (0) analytically. However, by change of variables from 
rj to a, these equations become remarkably tractable: 



E{a, a) 



''I drj 



da 



fi a{riY Ja a'^Va + 



(B2) 



and 



H(a, a) 



daa^ 



a y/a + a^ 



E(a, a) 



-(a — a) 



+ — vTT 
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n-3 



-5/6 



^4' 6' 6 ' 



(B3) 



Here, 2-^1 is Gauss' hypergeometric function. Gradshteyn and Rhyzik equation (3.194.1) was 
used to compute (BBf). Equations ([B2D and (|B^) can be used to compute GBR fluctuations 
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in a vacuum dominated model. Note that ( p2D and (p3| ) can also be used to evolve vector 
and decaying scalar modes: 



hl(n) = 16nG J Q^{f])E{a,a)a^df] 



APPENDIX C: CLOSED MODELS 



In a closed matter-dominated FRW model, the scale factor is 



a{ri) = 1 — cos(?7) 



(CI) 



Thus, the scalar modes are given by equation (|^) with 



6a 

H{v,f]) = (Y-2sin(r7) + 



cos(?7) sin(?7) , 



(C2) 



The vector modes still satisfy equation (|T3D and with the appropriate form for a{ri), equation 



(111) describes the evolution of tensor modes. 
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FIGURES 
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Poincare to Friedman mapping 

FIG. 1. The solid lines are the equispaced surfaces of constant z. The dotted lines are 
geodesies of constant x. The periodicity boundary is selected along one of these curves. The 
crosses indicates the location of our numerical grid points, which are regular and evenly spaced in 
the Poincare metric. 
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Photon Trajectories 

FIG. 2. This figure depicts various photon geodesies in the Poincare frame. The axes are in 
units of constant distance, so p = exp(— 2;)(a;^ + y^) and z = ln{w). The heavy fine is the fudicial 
geodesic along which our grid is periodic. The radial lines is the trajectory traversed by photons. 
The concentric lines are circles of constant distance from the origin, i.e. spheres in the Poincare 
frame. 
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Q = 0.21 Multipole Spectrum 

FIG. 3. The eight jagged lines correspond to observers in different universes or at different loca- 
tions. Flat inflationary models predict + 1) as constant. The dashed line is a parametrization 
of the CBR open universe spectrum: cil{l + 1) oc /(I + with ^max = 8 and q = 2.5. In 

general Zmax scales as oc 
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Q = OA Multipole Spectrum 

FIG. 4. Same as figure |, but for n = 0.4. 
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Q — 1 Multipole Spectrum 

FIG. 5. Nearly flat universe calculated using the open universe code. Here we recover the flat 
Harrison Zeldovich spectrum. 
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Flat Space Approximation 

FIG. 6. Flat space simulation at z=l 

FIG. 7. Flat space simulation at z=1.5 
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Comparison with APM survey 

FIG. 8. Comparison of the predicted open universe power spectrum to the APM survey for 
various parameters of Hubble's constant and CI 
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